Stability of the lattice formed in first-order phase transitions to matter containing 

strangeness in protoneutron stars. 
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Well into the deleptonization phase of a core collapse supernova, a first-order phase transition to 
matter with macroscopic strangeness content is assumed to occur and lead to a structured lattice 
defined by negatively charged strange droplets. The lattice is shown to crystallize for expected 
droplet charges and separations at temperatures typically obtained during the protoneutronstar 
evolution. The melting curve of the lattice for small spherical droplets is presented. The one- 
component plasma model proves to be an adequate description for the lattice in its solid phase with 
deformation modes freezing out around the melting temperature. The mechanical stability against 
shear stresses is such that velocities predicted for convective phenomena and differential rotation 
during the Kelvin-Helmholtz cooling phase might prevent the crystallization of the phase transition 
lattice. A solid lattice might be fractured by transient convection, which could result in anisotropic 
neutrino transport. The melting curve of the lattice is relevant for the mechanical evolution of the 
protoneutronstar and therefore should be included in future hydrodynamics simulations. 

PACS numbers: 26.50.+x,64.60.Cn,97.60.Bw 



I. INTRODUCTION 

In the past years, considerable effort has gone into the 
study of the existence of phases of strange matter in neu- 
tron stars jl], ^, ||, ^, g). A macroscopic strangeness 
content is predicted to be energetically favorable at den- 
sities well beyond the saturation density of symmetric nu- 
clear matter (hereafter p ). Conditions with sufficiently 
high density and low enough electron chemical poten- 
tial for the formation of a macroscopic content of strange 
quarks (charge = —1/3) or particles containing such are 
obtained in the interior of a protoneutronstar (hereafter 
PNS) after post-bounce times of several seconds. The 
three possible forms for the macroscopic manifestation 
of strangeness suggested are a K~ condensate (|, ||] , the 
formation of hyperons and deconfined quark matter 
including strange quarks [}) . The equation of state of 
matter at densities p > po is not well enough known at 
this point to determine which of these scenarios will ac- 
tually happen and what the order of the associated phase 
transition is. 

The present study limits itself to a certain class of 
scenarios which is common to all possible forms of 
strangeness in high density nuclear matter and which 
might be subject to experimental verification. I assume 
the phase transition to be first order and to result in 
the formation of a lattice in the coexistence region with 
different energy and charge densities of the strange and 
non-strange phases, as was recently predicted by various 
authors @, |, |, |g, 0, ||. Typical results for the pa- 
rameters of the mixed phase lattice, such as the size and 
spatial separation and the energy and charge densities 
are quoted from these authors. Whereas all studies on 



this lattice to date are limited to cold, deleptonized neu- 
tron stars, my focus is on PNS's in the Kelvin-Helmholtz 
cooling phase, ~ 1 s — 30 s post-bounce, for tempera- 
tures up to tens of MeV. In particular, the present pa- 
per presents the mechanical properties of the structured 
mixed strange/non-strange phase as they are relevant to 
the further PNS evolution and observational verification. 

Many initial studies predicted phenomena caused by 
strangeness to be confined to a very small central region 
of the PNS Jl^, [l4| , which would most likely render their 
experimental verification impossible. However, Glenden- 
mng showed in 1992 § that the previous studies were 
too simplistic in their assumptions about the thermody- 
namic model of the PNS. In particular, there are two 
conserved charges in PNS matter, the baryon number 
and the electric charge and, therefore, two independent 
chemical potentials, p e and ps- During a phase transi- 
tion, their equality between the strange (subscript s) and 
non-strange (subscript n) phases, together with similar 
conditions for temperature and pressure, comprise the 
Gibbs conditions j7]: 



PB,. 
P, = 



Pn. 
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Previously, charge neutrality had been enforced in both 
phases of a first-order phase transition separately, and a 
Maxwell construct had been used for the mixed phase: 
Pmixed = XPs + (1 - x)pn (x being the volume fraction 
of the strange phase) , where the densities of each phase 
(p s , p n ) were kept constant throughout the phase transi- 
tion. This additional constraint lead to an overestimate 
of the required density for a phase transition to decon- 
fined (three flavor) quark matter of up to 7 x po as op- 
posed to a predicted range of ~ 2 — 3 x p [0, || with 
the correct treatment of the Gibbs condition. The mi- 
croscopic stability condition for neutron stars SP/ dp > 
]15[ and the constant densities in both phases also lead 
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to a coexistence region consisting of an infinitely thin 
spherical shell with a pressure discontinuity, as opposed 
to a spatially extended mixed phase with continous P(r)- 
dependence which follows from the Gibbs conditions [[10| . 

The driving force for a net exchange of charge between 
the strange and non-strange phases is the isospin restor- 
ing force which leaves the strange phase with a net neg- 
ative charge. Two phases with opposite charges assume 
a spatial order which is determined by a minimal sum 
of the surface, curvature and Coulomb energies & [To|| . 
The predicted geometry of the strange phase has been 
found Jl^, [ll| to vary from spheres via rods to platelets 
immersed in the majority non-strange phase, as density 
increases. As soon as x > 0.5, the two phases reverse 
roles. The spatial extent of the crystalline phases de- 
pends sensitively on the neutron star mass [jl0| and prop- 
erties of matter at supernuclear densities which are not 
well known, but can reach several km. Some of these 
properties have been investigated in recent works, such 
as the surface energy between the strange and hadronic 
phases and the effective MIT bag constant in the case of 
a phase transition to deconfined quark matter ]T^ | and 
the surface and curvature energies between normal nu- 
clear matter and a phase with a K ~ -condensate |Q. In 
both cases, the results were shown to be mo del- dependent 
without a definitely reliable result. 

The above quoted range of ~ 2 — 3 x p |], || for 
the transition density is based upon the bulk approx- 
imation for both phases, which neglects any screening 
effects across the interface. More recent studies of sur- 
face effects, such as a transition layer with a finite thick- 
ness of AR ~ 5 fm due to Debye screening effects and 
the discontinous pressure due to the surface tension be- 
tween both phases, however, show a net increase in the 
bulk energy density of the strange phase in the case of a 
K~ condensate |l7j , making finite-size droplets for radii 
smaller than R$ ~ 10 fm energetically less favorable. De- 
bye screening lengths have been reported for K~ con- 
densed matter (A^- <~ 5fm [[171), for deconfined quark 
matter (X q ~ 5fm j|) and for hadronic neutron star 
matter (Xd,u/p ~ 10 fm and A e - ~ 13 fm ||). Properly 
taking into account screening to find the minimum en- 
ergy configuration, see for K~ condensates and Q 
for deconfined quark matter, effectively opens up another 
degree of freedom, in itself lowering the energies of both 
phases, which is favorable for the formation of a lattice. 
However, screening can also increase the electron frac- 
tion in the hadronic phase by pushing the electrons away 
from the negative charge on the strange droplet surface, 
leading to a higher negative charge concentration out- 
side the screened region which will push the pressure at 
which glob al charge neutrality can be attained to higher 
values |fl7f . The screened charges have also been shown 
to increase the effective surface tension a — <J s trong + °c 
by a Coulomb contribution, making droplets with radii 
smaller than the Debye screening length energetically less 
favorable (see Jl8| for deconfined quark matter), thus 
increasing the required PNS density for the first-order 



phase transition. A final answer is not possible unless 
the exact equation of state for hadronic PNS matter with 
and without a K~ condensate and deconfined quark mat- 
ter with strange quarks is known. This includes hyperon 
formation, following which a similar structured phase has 
not been studied to date. 

Standard PNS models predict central densities in the 
relevant range of a few times po after neutrinos have car- 
ried away the bulk of the lepton number in the PNS. 
A mixed strange/non-strange phase which might form 
then was shown to have a significant impact on the neu- 
trino transport properties, resulting in a neutrino opacity 
up to two orders of magnitude higher for typical neu- 
trino energies ~ lOMeV [flSfl . This opens up a window 
of observation in supernova neutrino detectors through 
the remaining lepton number which will be carried away 
by neutrinos. A pure strange phase which might form for 
high enough central densities would have a lower neutrino 
opacity compared to non-strange PNS matter p(| , effec- 
tively creating a transparent strange core surrounded by 
a relatively opaque coexistence layer. 

In the following, I will refer to a structured mixed layer 
as phase transition lattice, independent upon the specific 
spatial ordering and solid or liquid state. The temper- 
ature behavior of the phase transition lattice, including 
the coexistence curve between a solid, crystalline phase 
transition lattice and a liquid phase of the droplets con- 
taining strangeness is derived in section |l[ The mechan- 
ical stability of the lattice to shear stresses and its pos- 
sible breakup is investigated in section |l|. Section IV 
concludes the study and gives an outlook on some possi- 
ble observational signatures of the formation of the phase 
transition lattice. 



II. THE PHASE TRANSITION LATTICE AT 
FINITE TEMPERATURES 

A. Phase Transition Lattice Melting Curve 

Because a solid, crystalline phase might have to be 
taken into account as a new element in hydrodynamical 
PNS evolution studies and since its neutrino transport 
properties might be different from a liquid mixed phase, 
it is important to know the solid-liquid coexistence curve 
of the phase transition lattice. At the onset of the first- 
order phase transition and on the outermost layers of the 
mixed phase in the cold neutron star, the volume frac- 
tion of the strange phase will be small. In that limit, the 
minority phase at zero temperature can be approximated 
as a Coulomb lattice of negative point charges immersed 
in a slightly positively charged background of normal 
PNS matter. The absolute charge density in the strange 
phase will be higher because the majority phase can sig- 
nificantly lower its isospin by pushing negative charge 
into the minority phase p^ , more than compensating 
for the opposing effect of a higher repulsive Coulomb en- 
ergy within the latter. The reverse is true for the deepest 
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layers of the mixed phase, where the normal PNS matter 
is the minority phase and a high positive charge density 
resides in that phase, whereas the Coulomb interaction 
is compensated by the condensation energy of hadrons. 

The droplet charges given are to be understood as ef- 
fective values, already taking into account screening ef- 
fects H [I?], . Given rigid strange droplets with charge 
densities large compared to the surrounding hadronic 
PNS matter, the lattice can then be regarded as a one- 
component plasma (OCP). Its treatment in the harmonic 
approximation is well-established |2^, ^3], ^4], |25|. The 
equation of motion of the a-componcnt of the displace- 
ment u on a lattice site 1 for a Bravais lattice (defined as 
having one particle per unit cell) can be written as 



Mu a (l) 



5$ 

Su a (l) 



■E 

PV 



<5 2 $ 



8u a (l)5up{l') 



7TMO, (2) 



where $ is the electrostatic potential and the displace- 
ment amplitude u a is 



2NM 



i(k) 



(3) 



with Ag. 



ag. and the usual definitions for 

creation- and annihilation operators. The three charac- 
teristic polarization modes, two transverse and one lon- 
gitudinal, of the Bravais lattice are denoted by the index 
3- 

In this formalism, the mean square displacement rel- 
ative to the distance between nearest neighbors d = 
(37r 2 ) 1 / 6 a for a BCC (body-centered cubic) lattice is 



1 h ^coth(fl%u}j(k)) 
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(4) 



where rj is the degeneracy parameter and the Debye in- 
tegral is defined as 



D n (x) = — / dt{-—). 
x n Jo e* - 1 



(5) 



The dispersion relation for the transverse modes has the 
acoustic Debye form [pil p5l 



u>T(k) = CtLUp- 



(6) 



the average square displacement amplitude. For typical 
lattice constants of a ss 10 fm, we obtain kjj « 0.4 fm" 1 
and huj p w 5.8 MeV. Typical values for the strange 
(minority) phase were used here, a charge density of 
pc = 0.4 fm -3 , a mass density of pM = 0.4 fm -3 and 
a droplet radius of R = 3.0 fm. The plasma frequency is 
an important quantity characterizing the lattice, because 
the degeneracy parameter rj = hup/ksT determines the 
role quantum effects play and, ultimately, the freeze-out 
of the OCP. In the present case, for "typical" protoneu- 
tronstar evolution temperatures, T ~ 10 MeV, we get for 
the degeneracy parameter rj — (hZe)/(^MeokBTa 3 / 2 ) ~ 
0.5. The problem at hand can therefore be treated nei- 
ther in the zero temperature- (quantum-) nor in the clas- 
sical limit. 

The Lindemann parameter j 2 is defined as the value of 
the quantity (u 2 )/d 2 at the solid-liquid transition. For an 
OCP, it has been determined using Monte Carlo - sim- 
ulations in both the classical (high temperature) limit 
and in the quantum case (zero temperature) 
for both fermions and bosons. Energetically, spin 
pairing effects on the droplets will drive their total spin 
to zero. I therefore treat them as bosons. For the inter- 
mediate degeneracies given here, the interpolation for- 
mula for the Lindemann parameter of a bosonic OCP by 
Chabrier |25fl is used: 




= 7o ~ 



0.096 + 4.31 x 10" V 



1 + 0.05?/ 2 



2.092 x 10" V : 



(7) 



with 70 = 0.249 being the quantum limit. The melt- 
ing curves for different charge densities on the strange 
droplets are plotted in figure [I], where the curves rep- 
resent charge densities from pc = 0.1 fm~ 3 to pc = 
0.8 fm~ 3 in steps of pc = 0.1 fm~ 3 for a mass density 
of pm = 0.4 fm~ 3 and droplet radii R = 3.0 fm. If the 
droplets are treated as fermions, the result only differs 
significantly for coexistence curves with transition tem- 
peratures below 1 MeV. It can be seen that an initially 
liquid phase transition lattice crystallizes for PNS tem- 
peratures of T ~ 1 — 10 MeV, which lies well within 
the range obtained during the Kelvin-Hclmholtz cool- 
ing timescales in the standard PNS paradigm, for a wide 
range of parameters (% and pc)- There is, for any given 
charge density, a lower limit on the droplet radius be- 
low which there is no crystallization. For example for 
pc = 0.4fm~ 3 , this limit is between 1.5 — 1.75 fm, see 
figure [| Hence, lattices with larger droplet radii solidify 
earlier in the PNS evolution, causing the mixed layer to 
freeze out from its interior outwards. 



with the Debye wavenumber ko = (6n 2 N /V) 1 ^ , the 
plasma frequency uj p = (Ze/eo x N/MV) 1 ^ 2 and a = 
0.393. For the longitudinal branch, the Einstein model 
has been suggested with a constant frequency of ujl oc lo p 
p4j . However, for a cubic lattice in the harmonic ap- 
proximation, the symmetry condition (u 2 ) = 3 x {u\) = 
3 x {u 2 L ) |22| makes only one branch necessary to calculate 



B. Deformation Modes 

The OCP assumes negative, inherently rigid, point 
charges in a sea of positive background. However, since 
the surface tension is only a ~ 10MeV/fm 2 , small 
compared to strong interaction energy scales e s trong ~ 
10 3 MeV/fm 3 for the given droplet dimensions, it is clear 
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FIG. 1: Melting curve for different droplet charge densities: 
pc = 0.1 fm~ 3 (bottom curve) to pc = 0.8 fm -3 (top curve) 
with p = 0.4 fm -3 and R = 3.0 fm. 
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FIG. 2: Melting curve for different droplet radii with pc = 
0.4 fm" 3 and p = 3.0 fm" 3 . 



which, when expanded to second order in dR, yields 

S ~ AttR 2 + 2^dR 2 = S + AS, (10) 
8 

from which the elastic constant ks for the deformation 
energy can be deduced: 



AE S = aAS = 2ir-adR 2 = -k s dR 2 . (11) 



The inertial term ms can be found via the kinetic energy 
i-R r VR 2 -x' 2 



dx 
-R Jo 



dp{2*p[\p M x 2 (^-) 2 ]) 



= ^PM^ S R 5 = i(| PM R 3 )R 2 = \m 3 R 2 . (12) 

The characteristic vibration energy can therefore be es- 
timated as 



L0 S 




(13) 



A typical value, for pm — 0.4 fm~ 3 , R — 2fm and 
(7 = 10McVfm~ 2 , is hojs ~ lOMeV, which is compara- 
ble to the plasma frequency. Hence, deformation modes 
freeze out at about the same temperature as lattice vibra- 
tions. The OCP can therefore be considered as a valid 
description of the melting curve of the phase transition 
lattice, since no other modes are relevant once it becomes 
a crystal. At temperatures above the transition between 
a liquid and a solid strange droplet phase, the lattice- 
and deformation modes will be in thermal equilibrium. 



that the strange phase droplets cannot necessarily be 
considered as rigid. It is therefore important to know 
whether deformation modes have to be taken into con- 
sideration in the treatment of lattice vibrations. 

Consider a droplet of strange matter which is slightly 
elongated along the x-direction to R + dR, yielding an 
ellipsoid: 



Its surface area is 
S = 2n(b 2 



Va 2 - b 2 



■ /« 
arcsmy- 



-)), 



(8) 



(9) 



III. SHEAR STRESSES AND THE STRANGE 
PHASE TRANSITION LATTICE 

Besides thermodynamic criteria for the existence of a 
crystalline mixed phase, we need to know whether the 
lattice can withstand typical shear stresses present in 
convective PNS cores. The shear constant of a cubic 
Coulomb lattice is [B0| 



C44 



d 2 Wi 

dixy 



(14) 



where Wi is the lattice energy and ~f xy the angle of dis- 
tortion. The Coulomb lattice energy can be calculated 
using Ewald's method pi], p2] : 
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2 ^ 47re r(0 
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crfc(<?r(7)) 
r(0 



E 



4tt exp(-Gf/4g 2 ) 



G 2 



(15) 



r 



where the complementary error function erfc(x) = 
^= J x exp(—y 2 )dy, f2 is a unit cell volume and g is a 
parameter to be adjusted for fast numerical convergence 
of both real (vectors r(l)) and inverse (vectors G(l)) lat- 
tice sums. The result for a bcc-lattice is [BOl 



C44 



0.7423 x 



\7TR 3 p c e 2 
Aireoa 



= 4.477 MeV x (R 3 pc)( 



1.0 fm 



(16) 



a being the lattice constant. The critical shear stress is 
the force per unit area necessary to maintain two planes 
of the crystal distorted against each other by an angle 
corresponding to a displacement of a/4 perpendicular to 
a lattice plane |33[, which is in the linear approximation 



1 dll Id. r x \2\ 
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which, for a set of typical values, a lattice constant 
a = 10 fm, droplet radius R — 2.0 fm and charge den- 
sity pc = 0.4 fm~ 3 , gives a cnt = 1.4 x 10~ 3 McVfm~ 3 . 
Stresses in that order of magnitude might be caused by 
convection or differential rotation of the newly formed 
PNS. These stresses can either prevent the formation of 
the solid lattice in the first place or, if they are due to 
phenomena which are prone to variations, such as con- 
vection, break up a solid lattice formed during a transient 
period of weak convection. 

A negative gradient in the lepton concentration has 
been shown to lead to convection during the Kelvin- 
Helmholtz cooling phase of the PNS (3^] . More recently, 
hydrodynamics simulations including convection indicate 
that the Ledoux criterion for convective instability 



Cr 



Sp 

ss 



P.Y, 



dS_ 

dr 



dY, 



Sp 



SYi J p s dr 



> 



(18) 



is true in most of the PNS for times of more than ~ 1 s 
after bounce J35|. The pure strange phase in the center 
of the PNS, if it exists, is not expected to show strong 
negative lepton or entropy gradients. This is due to the 
relatively opaque mixed phase enclosing it and the fact 
that the transport of both heat and lepton number in 
the PNS interior is dominated by neutrinos. The most 
violent convection will therefore take place in the matter 
exterior to the phase transition lattice. Although various 
authors disagree on the extent and strength of convec- 
tion in PNS's, convective velocities of v c ~ 10 6 ms _1 are 



reported in many studies 35, jJTJ. This is equivalent 
to a kinetic energy density of E CO n^/V ~ 10~ 3 MeV/fm 3 , 
which indicates that convection might indeed be able to 
either break an existing phase transition lattice or pre- 
vent its formation. A strong enough convective cell form- 
ing outside the mixed strange/non-strange layer after 
its crystallization during a transient quiet period might 
fracture the solid lattice and mix matter from the non- 
strange envelope into the now liquid phase transition lat- 
tice in the region below the convective cell. For the dura- 
tion of the convective flow, this would result in a localized 
hole with a substantially lowered neutrino opacity JTs) ] 
compared to the still intact solid lattice in all other di- 
rections, hence in anisotropic neutrino transport through 
the mixed strange/non-strange layer. 

The discovery of a number of millisecond pulsars in 
recent years |3q ] indicates that some of the angular mo- 
mentum residing in the core collapse supernova might 
remain in the PNS. The resulting rotation is likely to 
be differential and has been studied by Goussard et al. 
p9| , |i"o| who solved the relativistic stellar structure equa- 
tions for rotating PNS's with representative equations of 
state for the different epochs in the PNS evolution. The 
rotation period £1 as a function of radius r assumed in 
that study is (in the Newtonian limit) Efl] 



n = 



R 2 £lc 



R 2 + r 2 sin 2 (9) 7 



(19) 



with Ro ~ 1 km the characteristic scale of variation of 
f2, r sin(8) the distance from the rotation axis and ilc 
the central rotation period. The rotation period a PNS 
can acquire without additional accretion has been shown 
to be limited by the increase of the minimum neutron 
star mass for times up to ~ 100 ms post-bounce and 
by the mass shedding limit beyond that, resulting in 
Pmin ~ 1.7 ms |4(J. This corresponds to velocities of 
v ~ (r/1 km)(ri/rTc , )10 6 ms _1 , which is comparable in 
magnitude to convective velocities sufficient to cause crit- 
ical shear stresses. It is not likely that regions of a fast 
rotating PNS at several seconds post-bounce go through 
a transient phase with low rotation period during which 
the phase transition lattice could crystallize, possibly to 
be broken up at later times by a larger gradient in the 
rotation period. Rather, for PNS's with rotation periods 
below ~ 100 ms, the timescale for the stratification of dif- 
ferential rotation might be of significant importance for 
the melting curve of the phase transition lattice, possibly 
comparable to the cooling timescale. However, unless the 
transport of angular momentum in core collapse super- 
novae and in particular within the PNS is finally resolved, 
the melting curves presented in sec. II A only apply to 
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PNS's with low differential rotation. 



IV. SUMMARY AND DISCUSSION 

The present study shows that if the phase transition 
from hadronic PNS matter to strange matter is of first 
order and if it results in a lattice of separate strange 
and non-strange phases in the coexistence zone, it will 
crystallize for temperatures predicted during the Kelvin- 
Helmholtz cooling phase of a PNS following a core col- 
lapse supernova. The process of crystallization has, how- 
ever, complex interactions with the hydrodynamical evo- 
lution of the PNS. Shear stresses due to strong differ- 
ential rotation with minimum periods in the range ob- 
served for millisecond pulsars and convection induced by 
lepton gradients might, as long as they persist, prevent 
the formation of a solid lattice. A final answer can only 
be given by the full inclusion of a possibly solid mixed 
strange/non-strange phase in future core collapse super- 
nova simulations. 

If the first-order phase transition occurs before the 
end of the Kelvin-Helmholtz cooling phase, the remain- 
ing neutrinos to be emitted might serve as a possible 
window on its formation and crystallization. With the 
predicted increase in the neutrino opacity for a first or- 
der phase transition lattice JnJ, expected to be espe- 
cially pronounced for intermediate neutrino energies in 
the range ~ 10 — 100 MeV, the timescale for the decay 
of the neutrino emission will increase, which might be 
observable as a knee in the neutrino luminosity. Once 
the lattice becomes solid, the evolution might be accom- 
panied by fractures and rearrangements in the mixed 
strange/non-strange zone, possibly showing irregularities 
in the neutrino luminosity. If a changing lepton number 
gradient causes a convective cell to form in a formerly 
non-convective region located outside a previously crys- 



tallized phase transition lattice, sufficiently violent con- 
vection might break the solid lattice and locally cause 
it to return to a liquid state. Non-strange PNS matter 
transported from regions outside the phase transition lat- 
tice might mix with the liquid, resulting in an anisotropic 
neutrino-opaque layer with a partial hole and therefore 
anisotropic neutrino transport. This will be more pro- 
nounced if the central density of the PNS is high enough 
to allow for a pure strange core enclosed by the relatively 
opaque mixed zone which essentially dams up neutrinos 
behind it. The extent to which the phase transition lat- 
tice will affect the neutrino luminosity and emission spec- 
trum will be the topic of a detailed transport study with 
a special focus on anisotropic neutrino transport in con- 
junction with a hydrodynamical treatment of the non- 
strange layers above the coexistence zone [[II]. This will 
also include the different neutrino transport properties of 
a solid versus a liquid lattice. 
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